GUIDANCE AND CONTROL OF TACTICAL MISSILES 






Thomas Alan Grote 









NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

GUIDANCE AND CONTROL OF TACTICAL MISSILES 

by 

Thomas Alan Grote 
December 1979 

Thesis Advisors D. J . Collins 

Approved for public release? distribution unlimited 



T1961.76 



security classification of this mgk rwhm * Data em«r«dj 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


rerort number 


2. GOVT ACCESSION NO* 


l. REClPltNT’l CATALOG NUM8CA 


4 TITLE (and Subtitla) 

Guidance and Control of Tactical Missiles 


»• type or *tp o<*r * *e<*ioo covered 

Master’s Thesis; 

(December 1979) 


S. RERFORMING D RG . RERORT NUMBER 


7. author^ 

, Thomas Alan Grote 


s. contract or grant number^ 


• PERFORMING ORGANIZATION NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 939^0 


io rrogram element project, task 

AREA A WORK UNIT NUMBERS 


II CONTROLLING OFFICE NAME ANO ADDRESS 

Naval Postgraduate School 
Monterey, California 939^0 


12 REPORT DATE 

December 1979 


11 Nl^fcJ^ER OF P AGES 


T4 MONITORING AGEnCy name 6 AOO»ESSff/ ditfarant from Controlling Qtttca) 

Naval Postgraduate School 
Monterey, California 939^0 


IS. SECURITY CLASS for rfti# rm&ott) 

Unclassified 



Tia, DECLASS|F|CATIOn’downGRAO;nG 
SCHEDULE 



16 DISTRIBUTION STATEMENT ( ot thl a Ro^orf) 

Approved for public release; distribution unlimited. 



17. DISTRIBUTION STATEMENT (ot r*« amatrad antarad in Block 30, It dlffarant from Rapori) 



IS. SUPPLEMENTARY NOTES 



IS. KEY WORDS (Conltnum on ravaraa old* If nacaaamry and identity by block numbmr) 

M0D6DF - STM 



20. ABSTRACT (Confirm* on ravaraa aid* It nacaaamry and ld*nflty by black m^b*r) 

This thesis discusses the conversion of the M0D6DF computer program for use on 
the I3M-360 computer at the Naval Postgraduate School* The functioning program 
was modified to investigate the impact miss distance for the Supersonic 
Tactical Missile. When the initial y-displacement error exceeded 1800 feet, 
the missile did not acquire the target. All errors smaller than this resulted 
in miss distances within 0.5 feet of the target. The midcourse guidance 
reference altitude was changed to reflect a sea -skimming missile. This 



DD 1473 

(Page 1) 



EDITION OF t MOV «» IS OBSOLETE 
S/N 0 103-0 14- 660 1 j 



1 



SECURITY CL ASSI FICATlON OF THIS RACE (Whan D*t* Bntarad) 



Cl*>»i>iC*^iON Q» nrnim Mmtmr+a 



simulation ran and impact was recorded. An attempt at adding random noise 
to the homing seeker was tried, but revealed that more information is required 
on this topic. The M0D6DF computer program was successfully converted and al- 
tered to run using the simplified ramjet model. 



DD 

S/ 



Form 1473 
1 Jan ,3 

N 0102-014-6601 



2 



ticuxvv classification of tn»i ^atr**** a*'* *-'*'**> 



Approved for public release; distribution unlimited 
Guidance and Control of Tactical Missiles 

by 

Thomas Alan Grote 
Lieutenant, United States Navy 
B.S.A.E., United States Naval Academy, 197*+ 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER UF SCIENCE IN ENGINEERING SCIENCE 



from 



NAVAL POSTGRADUATE SCHOOL 
December 1979 



Crt! 0>C 



ABSTRACT 



This thesis discusses the conversion of the MCC6DF computer program 
for use on the IBM- 36 O computer at the Naval Postgraduate School. The 
functioning program was modified to investigate the impact miss distance 
for the Supersonic Tactical Missile. When the initial y-displacement 
error exceeded 1800 feet, the missile did not acquire the target. All 
errors smaller than this resulted in miss distances within 0.5 feet of 
the target. The midcourse guidance reference altitude was changed to 
reflect a sea -skimming missile. This simulation ran and impact was 
recorded. An attempt at adding random noise to the homing seeker was 
tried, but revealed that more information is required on this topic. 

The M0D6DF computer program was successfully converted and altered to 
run using the simplified ramjet model. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 



The following is a list of the abbreviations and some of the more 
common fortran symbols used in the M0D6DF computer program. Each symbol 
is defined in two ways. The primary source of identification is by its 
COMMON(3M5) location and the second is by its fortran symbol. The 
fortran symbol is not always a good identifier since it will change from 
subroutine to subroutine (i.e. W and WE both are used for missile 
weight) . 



ABBREVIATIONS 










Symbol 


Definition 








TP 


tangent plane axes 








BA 


body axes 








SA 


stability axes 








FORTRAN SYMBOLS 










Symbol 


Definition 








TXBA(073) 


missile thrust in the x-direction 


in 


BA 


(lb) 


TYBA(07^) 


missile thrust in the y-direction 


in 


BA 


(lb) 


TZBA(075) 


missile thrust in the z-direction 


in 


BA 


(lb) 


W( 086) 


missile weight (lb) 








S(110) 


missile reference area (ft ) 








CA(lll) 


drag coefficient 








CY(113) 


side-force coefficient 








CZ(115) 


normal force coeffieient 








CBAR(ll6) 


mean aerodynamic chord (ft) 








CMQ(ll8) 


damping in pitch coeffieient 








CNR(H9) 


damping in yaw coefficient 
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Symbol 

CLP(l20) 

CM(l2l) 

CN(122) 

CL(123) 

CGI(136) 

A(20l) 

B(202) 

CC(203) 

TSA(208) 

P(212) 

Q(2l6) 

R(220) 

AG (282) 

VXTP(286) 

XTP(290) 

VYTP(29^) 

YTP(298) 

VZTP(302) 

ZTP(306) 

GZRO(404) 

ER(405) 

ALPO(406) 



Definition 

damping in roll coefficient 
pitching moment coefficient 
yawing moment coefficient 
rolling moment coefficient 
center of gravity (ft) 

moment-of-inertia about missile roll axis (x-body 
axis) (slug-ft^) 

moment-of-inertia about missile pitch axis (y-body 
axis) (slug-ft^) 

moment-of-inertia about missile yaw axis (z-body 

axis) (slug-ft^) 

angle between SA and EA (rad) 

missile angular velocity about x-BA (rad/s) 

missile angular velocity about y-3A (rad/s) 

missile angular velocity about z-EA (rad/s) 

unit conversion lb-slug 

missile velocity in x-Tr (ft/s) 

missile displacement in x-TP (ft) 

missile velocity in y-TP (ft/s) 

missile displacement in y-TP (ft) 

missile velocity in z-TP (ft/s) 

missile displacement in z-TP (ft) 

constant set to zero for flat earth gravitational 
field and set to one for a spherical gravitational 
field 

angular velocity of earth (rad/s) 
angle between north and x-TP (rad) 
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Symbol 

0LAM0(407) 

H0(4l4) 

GO(4l5) 

HREF(50l) 

RE (503) 
H(507) 
AMACH( 520) 
THET(521) 
?SI(522) 

PHI(523) 

GAMMAV(527) 

VEL(528) 

VAT(529) 

TF(550) 

VAH(56l) 

AIAT(576) 

AZ(578) 

DYNP(58l) 

FLGRJ(606) 

T(932) 

Tl (933 ) 
T2(934) 
T3(935) 
T4(936) 
T5(937) 



Definition 

latitude origin of tangent plane 

distance tangent plane is from earth (ft) 

gravitational acceleration (ft/s^) 

reference altitude for midcourse guidance (ft) 

earth's radius (ft) 

altitude normal to earth (ft) 

missile Mach number 

missile pitch angle, TP (rad) 

missile yaw angle, TP (rad) 

missile roll angle, TP (rad) 

vertical flight path angle, TP (rad) 

magnitude of missile velocity (ft/s) 

missile velocity (ft/s) 

program termination time (s) 

computed speed of sound (ft/s) 

latitude of target position (deg) 

azimuth of target position (deg) 

dynamic pressure (psi) 

constant set to zero indicates run will use 
ENGINE subroutine and when set to one run will 
use RAMJET subroutine (simplified ramjet model) 

actual time (s) 

boost engine ignition (s) 

commence acceleration command mode (s) 

boost engine burn-out/port cover blow- in (s) 

start ramjet engine (cruise) (s) 

commence heading and altitude guidance control (s) 
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Symbol 

T6(938) 

T7(939) 

T8(940) 

DMAX(1740) 

CPP(2669) 

R0LL0(2901) 

PITCH0(2902) 

IAWO(2903) 

STEP (2905) 

DOC (2909) 

HMIN(2911) 

HHAX(2912) 

DER( 1)(2913) 



Definition 

commence terminal dive (s) 

commence terminal guidance - search (s) 

commence terminal guidance - track (s) 

maximum fin deflection (rad) 

time between printouts (s) 

initial roll angle (rad) 

initial pitch angle (rad) 

initial yaw angle (rad) 

determines executive program flow after staging 
defines number of times COMMON will be printed 
minimum integration step size 
maximum integration step size 
integration step size (s) 
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I. INTRODUCTION 



Research was undertaken to convert the M 0 D 6 DF computer program re- 
ceived from NWC China Lake for use on the IBM-36O computer at the Naval 
Postgraduate School. Once a functioning program was obtained, initial 
displacement error versus impact miss distance was investigated for the 
supersonic tactical missile. Additionally, the midcourse guidance 
reference altitude was changed to reflect a sea-skimming scenario and 
this was examined for its effect on the terminal guidance problem. 

This report not only discusses the aforementioned topics, but also de- 
scribes the missile mission requirements and the MGD6DF computer program. 



13 



II, MISSION REQUIREMENTS 



The Supersonic Tactical Missile (STM) mission is divided into six 
phases (Ref. l). 

* initial conditions 

* separation 

* boost 

* transition 

* cruise 

* terminal 

These divisions are based on the missile aerodynamics. The initial 
condition phase establishes the starting conditions for each launch. 

This is done while the missile is still attached to the launch platform. 
The separation phase starts when the missile is launched. The missile 
falls for approximately five seconds until the boost engine ignites. 

This initiates the boost phase which continues until the missile achieves 
Mach two. As the missile passes through Mach one, plume effects are 
encountered which the aerodynamics account for. At the end of th9 boost 
phase, the port covers blow in. This initiates the transition phase. 

This phase is very short and allows the debris to be ejected from inlet 
ports. The cruise phase commences when the ramjet engine ignites. This 
engine propels the missile until target impact. The terminal phase of 
the flight begins when the missile is commanded to dive from the cruise 
altitude. This phase concludes when the flight is terminated at target 
impact. 

The six STM mission phases require four phases of control. These 
control phases ares 



& 



* separation 

* midcourse guidance 

* terminal dive 

* terminal guidance 

The four Guidance, Navigation, and Control (GN<*C) system control phases 
are directly related to the mission phases* Figure 1 displays this 
relationship (Ref. 2). To obtain a successful flight, various types of 
control and guidance functions are required. ^Figure 2 (Ref. 3 ) illus- 
trates the guidance control mode sequencing in relation to the mission 
phases and also indicates the critical missile switching times. The 
control phases are discussed below, along with the appropriate guidance 
modes . 

A. SEPARATION 

The separation phase commences upon *launch. The missile is ejected 
downward from the launch platform. Additionally, the pitch attitude of 
the missile is commanded down. Since this portion of the flight is 
unguided, the ejection force and gravity are the only forces acting on 
the missile. At launch, the missile is required to be in the attitude 
command mode. 

Five seconds into the flight, the missile pitch attitude is commanded 
up and the boost engine is ignited. The booster continues until the 
missile attains Mach two, at which time control is shifted from attitude 
to acceleration control mode. The acceleration control then requires 
the missile to maintain the normal and lateral accelerations at zero. 

The last event to occur in the separation phase is inlet port cover 
blow-in. The time delay that allows the port covers to clear is a func- 
tion of altitude. 



15 



SEP BST TRN CRUISE TERMINAL 




16 



CONTROL PHASES /MISSION PHASES 
RELATIONSHIP 



to - LAUNCH: START SIMULATION 



ti - BuCSTER IGNITION 

t2 - CONTROL SYSTEM MODE CHANGE FROM ATTITUDE TO 
ACCELERATION CONDITION AT M=2.0 

tj - PORT COVER 3LCW-IN , BOOSTER BURNOUT 

t4 - RAMJET IGNITION, tj + .3 sec 

t5 - ENGAGE FUIDANCE MODES 

t6 - DIVE COMMAND 

ty - SEEKER SEARCH MODS IS ACTIVATED 
tg - TERMINAL TRACKING 



* 


IC 


SEP 


3ST 


TRN 


CRUISE 


TERMINAL 

| 1 — - 






t 0 i 


l i t' 2 




b ^9 ^3 

t 



CONTROL FUNCTIONS 

1. Attitude rate 

2. Acceleration 

GUIDANCE FUNCTIONS 

1 . Open loop 

2 . Heading 

3 . Altitude 

4 . Terminal dive 

5. Terminal guidance 

6. Terminal guidance 



search 

track 
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FLIGHT MODE SEQUENCING 
FIGURE 2 
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B. MIDCOURSE GUIDANCE 

Midcourse guidance begins upon completion of the separation phase. 

The moment the port covers are clear, the phases shift. During this 
phase two things are required. The GN&C must guide the missile to the 
established cruise altitude and then guide the missile (in the horizontal 
plane) to the predetermined target location. 

To accomplish these requirements the GNatC system employs both 
altitude and heading control. The altitude control acts upon the pitch 
axis to drive the missile to the cruise altitude and then to maintain 
that altitude until the terminal dive phase commences. 

The actual guidance work is performed by the yaw axis. The GNceC 
system uses the heading control to steer the missile to the target posi- 
tion. To generate the steering commands a guidance law is necessary. 

The guidance law should minimize the cross -track error and the final 
lateral acceleration. A minimum cross -track error allows for minimum 
flight time and minimum displacement error when target search is initi- 
ated. Minimizing the final lateral acceleration ensures that tne seeker 
will continuously be pointing toward the target area. To insure 
accuracy, the guidance law must be able to perform these functions when 
subjected to disturbances such as wind and thrust misalignment. 

C. TERMINAL DIVE 

The terminal dive phase commences when the missile is commanded to 
dive. The actual time that this command is given is a function of the 
predetermined target coordinates. Upon diving, the missile must accel- 
erate downward along a 60-degree (from horizontal) dive angle to the 
target position. Once the missile has steadied up in the proper dive 
aspect, terminal guidance commences. 
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D. TERMINAL GUIDANCE 

During the terminal guidance phase the seeker is required to search, 
acquire, and track the target. The GN<&C system then uses these tracking 
signals to direct the missile to the target. The seeker uses a prepro- 
gramed search pattern to locate the target. When the target has been 
acquired, the seeker then commences tracking it and also notifies the 
GN&C system that the target is being tracked. 

Target acquisition takes place when the target falls within the 
instantaneous four-degree beamwidth pattern as it traverses the scan 
pattern on the earth's surface. The computer program uses this assump- 
tion since the exact missile target acquisition mechanism has not yet 
been defined. 

Upon acquisition, the GN&C system closes the seeker tracking loop. 
The radiometer error signal is defined as the difference between the 
look vector and the line-of-sight . This signal is used to reposition 
the antenna gimbals so the look vector and the line-of-sight coincide. 
Under the closed loop operation the radiometer output is proportional 
to the line-of-sight rate and this is the signal used to guide the 
missile to target impact. 
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III. COMPUTER SIMULATION 



The modularized six-degree-of-freedom (MOD6DF) computer program :*as 
developed by the Litton Systems, Inc., Guidance and Control Systems 
Division (Ref. 4) to be used in analyzing missile guidance and control. 
The program uses a building block approach, where each module corresponds 
either to a missile subsystem or an environmental system. In its 
original form, the program is used primarily for terminal guidance of 
air-to-surface missiles. The Naval Weapons Center, China Lake modified 
the original program to specifically apply to the Supersonic Tactical 
Missile (STM). The modification also allows the user the capability of 
simulating any portion of the missile flight trajectory. 

The M0D6DF computer program consists of four main decks and one 
auxiliary deck. The main decks include the executive programs, 
operational subroutines, modules, and input data deck. All the sub- 
programs utilized in the integration and the sequencing of the modules 
are contained in the executive programs . The operational subroutines 
are used by the user to control the program while it is running. The 
psysical system and the environmental subroutines are included in the 
modules . The input data deck contains all the information the executive 
programs need to execute the desired run. Lastly, the auxiliary deck 
consists of all those subprograms (subroutines and functions) that are 
required by the modules . 

Once the user starts making simulations, he must concern himself 
with program sequencing. Proper sequencing is required to ensure a valid 
run is achieved. Since it is of such importance, sequencing will be dis- 
cussed as the final part of this section. 
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A. EXECUTIVE PROGRAMS 



The executive programs are the main core of the M0D6DF program. 

These subroutines have various functions in setting up, sequencing, 
integrating, and resetting the program. Since these functions are 
required by any type of analysis done, the user should never have to 
change any of these subroutines. Reference 5 should be consulted if the 
user desires more information about the content of any of the subroutines 
described herein. 

1 . Zero 

Subroutine ZERO is used to set all the COMMON ( 3 ^ 15 ) locations to 
the value zero. This is done to ensure that no erroneous information is 
used in the simulation. 

2. Oinptl 

Subroutine 0INFT1 is the basic input routine for the M0D6DF 
program. The normal input is from punched cards. However, inputs may 
also be read from tapes. The basic input cards will be discussed in the 
section covering the input data deck. 

3. Auxi 

Subroutine AUXI is used to call the initialization modules for 
each run. Additionally,, it sets up the list of state variables which 
are used in AMRK. 

4. Auxsub 

Subroutine AUXSUB is used to call the dynamic modules. In 
calling the dynamic modules, AUXSUB sets and resets the lists needed for 
AMRX and the COMMON(34l5) storage cells with the most recent values of 
the state variables and their derivatives . 



21 



5. Amrk 



Subroutine AMHK is the integration subroutine. It uses a point- 
wise first order Runge-Kutta method. All the state variables to be 
integrated must be listed in COMMON(34l5) and also must appear in a 
processing list. 

6. Reset 

Subroutine RESET is used to reinitialize up to fifty input 
parameters. This is done prior to the start of any repeated runs. To 
indicate which parameters are to be reset, the number one is punched in 
columns 46-60 on the respective type three card. 

7. Return Group 

The return group is a collection of all the unused modules . 

Since all the modules and their initialization modules are called by 
AUXI and AUXSUB, the unused ones must still remain in the deck. These 
are required to ensure proper linking when the computer attempts to 
link all the subroutines . All the subroutines in this grouping contain 
three cards. The three cards are the subroutine title, and a return and 
an end card. 

8. Subll, Subl2, Subl3 

These subroutines are used to call the operational subroutines 
that are required. They call the routines in the order prescribed by 
the input data deck. The number at the end of each subroutine title 
indicates which operational subroutines it can call. For example, 

SUB 11 can call STGE1 . 
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B. OPERATIONAL SUBROUTINES 

The operational subroutines provide the user with control of the 
program while it is running. The order in which the operational sub- 
routines are called is specified by the input data deck. Since these 
routines assist the user in controlling the simulation, they can be 
reprogrammed. However, it is advised that they not be changed until the 
user has become quite familiar with the overall operation and sequencing 
of the M0D6DF program. For more in depth knowledge of the operational 
subroutines the user should consult Ref. 6, and the computer listing, 
which is contained at the end of this report. 

1 . Inptl, Inpt2, Inpt3 

These subroutines are available for new inputs during the simu- 
lated flight. The only one presently used in the program is INFTl. It 
is utilized to input a namelist file, which is used by ENGINE and is 
described there. 

2. Oupt 1 , 0upt2 , 0uot3 

These subroutines allow for the print-out of up to fifty dif- 
ferent variables during the flight simulation. CUFTl is not utilized 
in the deck. 0UPT2 is used when the desired output is to be put on 
tape. 0UPT3 is the basic output routine. It prints the desired output 
on regular computer paper. 

3. Stgel, Stge2, Stge3 

These subroutines allow for proper staging, run termination, 
etc.. Presently STGSl is not being used. STGE2 is being used as the 
staging initialization subroutine. STGS3 is the primary subroutine of 
this group. It stages when impact with the earth is made, when the 

final time, TF(550)» is reached, or when LC0NV(2672) is set equal to 
two. All the tolerances for staging are listed in STGS3. 



23 



4. Cntrl, Cntr2, Cntr3 



These subroutines allow the external dynamic control inputs to 
the modules. In the M0D6DF program these are not used. 

5* Rndml, Rndm2, Rndm3 

These subroutines allow random noise to be added to the state 
variables generated in the modules. RNDM1 is not used in the program. 
RHDM2 is used as the initialization subroutine, while RNDM3 provides 
continous noise values. These subroutines produce correlated noise 
values for as many modules as required. The noise values remain fixed 
during each individual integration cycle. 

6. Auxal, Auxa2, Auxa3 
Auxbl, Auxb2, Auxb3 
Auxcl, Auxc2, Auxc3 

These subroutines are auxiliary routines that allow for external 
input, output, control, etc. of the modules. At the present time, none 
are utilized in the M0D6DF program. 

C . MODULES 

The modules are of prime importance to the user since they represent 
the f model f of the dynamic system. In general, the model is described 
by ordinary non-linear time-varying differential equations with both 
random and deterministic forcing functions. The user must first reduce 
these equations to an equivalent system of first order equations, which 
can then be described by each module. Generally speaking, the physical 
system is so complex that this would be impossible to do. However, due 
to the modularity of the M0D6DF program, the user can think of each 
module as a completely independent system described by the equations 
within that module. 

There are thirty-six possible modules divided into five functional 
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categories. Each group is identified by a letter which pertains to that 
groups function; A (airframe), C (computers), D (dynamics), G (geo- 
physical), S (sensors). A complete printing of each module is contained 
in the computer program listing. 

1 . Airframe 

a o Subroutine A1 

This is the aerodynamic forces and moments modules. It cal- 
culates all the necessary forces and moments in body axes. These 
values are then used in the computation of the dynamics. 

b. Subroutine A2 

This is the missile aerodynamic coefficient module. It cal- 
culates the required coefficients using the information stored in the 
BLOCK DATA. Using the timing inputs, this routine computes the coef- 
ficients for the different effects. Some of the effects accounted for 
are; plume, separation, and control surfaces effectiveness. With this 
done the total coefficients are determined. 

c. Subroutine A3 

This is the missile propulsion module. The timing inputs are 
used to determine whether the missile is in free fall, boost, transition, 
or cruise phase. With this determined the correct engine subroutine 
(BOOST, RAMJET, ENGINE) can be called. Three variables are calculated 
using the body axes as the frame of reference. They are the missile 
thrust is all coordinate directions, the principal moments of inertia, 
and the missile weight. 

d. Subroutine A4 

This is the fin actuator module. The four control surfaces 
commands are calculated as either ideal actuators or as second-order 
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ones. In addition to control surfaces commands, the rate of change of 
yaw, and roll are computed. 

e. Subroutine A5 

This module is part of the return group „ 

2 . Computers 

a* Subroutine Cl 

This is the autopilot module for the STV-G. It uses the 
cruise engine ignition time to divide the routine into boost /separation 
and cruise phases. These two phases use different algorithms to 
calculate the turning moments in pitch/yaw and roll. 

b. Subroutine C2 

This is the guidance command module. Using the timing inputs, 
this routine is divided into separation/boost, dive/climb, cruise, 
terminal dive, and terminal homing sections. Each section uses slightly 
different algorithms to calculate the guidance commands to maintain the 
proper flight profile. 

c. Subroutines C3 - CIO 

These modules are part of the return group, 

3. Dynamics 

a. Subroutine D1 

This is the translational dynamics module. It computes the 
total acceleration in body axes and then converts them to the tangent 
plane reference. Then, accounting for aerodynamics, thrust, gravity, 
and coriolis, the velocity and acceleration are calculated. 

b. Subroutine D2 

This is the rotational dynamics module. With the principal 

axes as a reference, this subroutine computes the body angular rates and 
the attitude direction cosines. 
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c. Subroutines D3 - D5 

These modules are part of the return group. 

4. Geophysical 

a. Subroutine G1 

This is the gravitational and coriolis acceleration module. 
It calculates the gravitational acceleration using one of two fields. 

The user specifies the field to be used by an input card. To use a 
flat-earth gravitational field, GZR0(404) must equal 0.0, and to use a 
spherical gravitational field, GZR0(404) must equal 1.0. 

b. Subroutine G2 

This module is part of the return group. 

c. Subroutine G3 

This is the air data module. It computes the velocity, in 
all three coordinate directions, with respect to the air mass. These 
values are then resolved into body and stability axes. This module also 
computes all the properties of air by calling subroutine AIR. These 
values are stored in their COHi , 'iON(34l5) locations for use in the other 
modules . 

d. Subroutine G4 

This module is part of the return group. 

e. Subroutine G5 

This is the coordinate conversion module. It takes the 
missile position, does a coordinate conversion and then it determines 
the position, velocity, and acceleration in the ECI system. 

f. Subroutine G6 

This module is part of the return group. 
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5. Sensors 



a. Subroutine Si 

This is the homing seeker module. It simulates the missile 
seeker and computes the seeker dynamics and £uler angle rates. Several 
flags are used to control the seeker sequencing: 

(1) FLAGS (33 5 ) 0 FLAGS signals the start of the seeker 

search. 

(2) FLGT ( 3 1 7 ) » FLGT signals when the seeker is locked-on. 

( 3 ) FLGTS( 3^7) <> FLGTS signals the end of search# 

(4) FLGD(336) q FLGD signals when the seeker has detected 

the target# 

( 5 ) FLAGLT . FLAGLT signals when the target is outside the 
seeker field of vision. 

b 0 Subroutine S2 

This is the radiometer module. It takes the target position 
and the ATIGS target position and converts them from body axes to the 
seeker axes. Using target position, the module then calculates the 
azimuth and elevation error signals. 

c. Subroutines S3, S4 

These modules are part of the return group. 

d. Subroutine S5 

This is the accelerometers and gyros modules. The user has 
the option of using Ideal or digital accelerometers. To specify the 
type of accelerometer, the user must include the appropiate data 
statement in the subroutine. If ideal accelerometers are desired, FLGA 
must equal 0.0 and for digital accelerometers, FLGA must equal 1.0. 

e. Subroutines S 6 - S10 

These modules are part of the return group. 
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D. INPUT DATA DECK 



The input data deck provides the user with the means of specifying 
which operational subroutines and modules are to be utilized for the 
desired run. It also allows the user to set the starting conditions. 

In general, only the constants and the state variables must be given 
initial values. All quantities in COMMON not given initial values will 
be set to zero by ZERO. In addition to the state variables, the upper 
and lower error bounds must be initialized. There are seven types of 
input cards, each indication a certain function. 

Type Function 

0 read/write tape 

1 operational subroutine to be called 

2 module to be called 

3 numerical input 

4 printed output 

5 parameter square and sum 

6 termination and random noise generator input 

A separate card is required for each subroutine, module, input, and 
output quantity. A sample computer printout of the input data deck is 
contained in the computer output section. 

1 . Type 0 Card 

Type 0 cards are used to indicate if the type 3 inputs are to be 
read from or written onto an auxiliary tape. These procedures can be 
used rather than reading the inputs from a deck of cards . A typical card 
is defined by punching a zero in column 2. The field that covers 
columns 5 “ 20 is used by the user for any descriptive statements with 
which he wishes to identify the input. Column 21 -23 contains the right- 
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justified integer number of the tape transport to be used. The number of 
the first record to be read is punched in column 31 * 45* The last 
field is column 46 - 60 which contains the number of records to be read. 
The last two fields may use either fixed or floating-point notation. 

2. Tyne 1 Card 

Type 1 cards are used to specify which operational subroutines 
are called during the flight simulation. This type card is identified 
by the number one in column 2. The second field is column 5 " 20 which 
contains any identifying information. This information is printed out 
when the data deck is read and allows the user to read exactly which 

subroutines were called. Column 21 - 25 contains the right- justified 
integer number which is the subroutine identifying number. The opera- 
tional subroutine numbers are; 

Subroutine Subroutine Number 

INPT1, INPT2, INPT3 2 

OUPTl , 0UPT2, 0UPT3 3 

STGEl, STGE2 , STGE3 4 

CNTR1, CNTR2, CNTR3 5 

RNDM1, RNCM2 , RNDM3 6 

AUXA1, AUXA2, AUXA3 7 

AUX31, AUXB2, AUXB3 8 

A0XC1, AUXC2 , AUXC3 9 

It should be noted that all or any of the subroutines listed under one 
one number can be called by including only one card. The cards are 
placed in the data deck in the order in which they will be called . 

This order or sequencing will be explained further in section F, 
Sequencing. A typical type 2 card is identified by the number two in 
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column 2* In general, column 5 * 20 should contain the module title, 
but any pertinent information is allowed. The module number is punched 



in column 21 - 25 and it 


must 


be right- justified . 


A listing of the 


module numbers follows : 








Module Module Number Module 


Module Number 


A1.A1I 


2 


D4.D4I 


20 


A2.A2I 


3 


D5.D5I 


21 


A3.A3I 


4 


• G1.G1I 


22 


A4,A4I 


5 


G2.G2I 


23 


A5.A5I 


6 


G3.G3I 


24 


C1.C1I 


7 


G4,G4I 


25 


C2.C2I 


8 


05,051 


26 


C3.C3I 


9 


G6 , G6l 


27 


C4.C4I 


10 


S1.S1I 


28 


C5.C5I 


11 


S2.S2I 


29 


C6.C6I 


12 


S3.S3I 


30 


C7.C7I 


13 


S4.S4I 


31 


C8 , C8I 


14 


S5,S5I 


32 


C9,C9I 


15 


S6,S6I 


33 


C10.C10I 


16 


S7.S7I 


34 


D1.D1I 


17 


S8,S8I 


35 


D2.D2I 


18 


S9.S9I 


36 


D3.D3I 


19 


S10.S10I 


37 


Note that either the module, 


the initialization module, or both may be 


called by including only 


one 


card in the deck. A 


sample of a typical 



type 2 card is shown in Figure 3 (Ref. 7). 
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COLUMN 
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5 20 


21 25 
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TYPICAL TYPE 2 CARD 
FIGURi 3 



2 




5 2C 


21 25 




31 ' 45 


46 60 




/ 3 




WEIGHT 


86 




1152.0 

1 .152 E+03 


1.0 

1.0 E+00 





TYPICAL TYPE 3 CARD 



FIGURE 4 



4. Type 3 Card 



Type 3 cards are used to set any C0MM0N(34l5) location to any 
value other than zero. In general, four items must be initialized. The 
state variable initial values and any constants used in the flight simu- 
lation are the most obvious. Additionally, there are some constants 
associated with the executive programs and operational subroutines and 
the state variable upper and lower bounds which must be initialized. As 
with all cards, column 2 defines the type card and it must contain a 
three. Column 3 * 20 holds the statement describing the input. The 
user should be specific here since it will save him having to remember 
every COMMON(34l5) location. The only other means of input identifica- 
tion is by column 21 - 25* These columns contain the right- justified 
COMMON location of the input. Columns 31 - 45 hold the actual numerical 
value of the input. The last field, column 46 - 60, contains the reset 
flag. If the reset flag equals one, the C0MM0N(34l5) location and the 
numerical value are placed in the reset list. This list may contain up 
to fifty different variables. This, in the case of multiple runs, 
allows the variables to be reset to its initial value prior to each run 
without additional input cards. The sample card in Figure 4 (Ref. 8) 
shows that either fixed or floating-point notation may be used to input 
the numerical value and the reset flag. These cards need not be inputed 
in any specific order, but for ease of checkout, it is advised to place 
them sequencially by their COMMON location. 

5* Type 4 Card 

The M0D6DF program can printout a maximum of fifty variables for 
each simulation. Type 4 cards are used to specify which variables are 
to be printed. Column 2 must contain the number four to indicate a type 



33 



k card. When the results are printed out headings aro included. These 
headings are designated in column 9 - 20. The exact alphanumeric title 
punched will be printed at the top of each page, this need not be the 
fortran symbol used within the program. The COMMCN(3M5) location of the 
output variable is contained in column 21-25 and must be right- 
justified. 

6 . Type 5 Card 

Type 5 cards are used to indicate which variables are to be root- 
mean-squared. These cards are similar to type 1 and type 2 cards. 

Column 2 must contain the number five. Any pertinent information about 
the variable to be operated on is punched in column 5 “ 20. The COMMON 
( 3 ^ 15 ) location of the variable must be right- justified in columns 21 - 
25* The last field, column 31 " ^ 5 , indicates whether the root-raean- 
square operation will occur along the trajectory or at the end. 

7. Type 6 Card 

The type 6 card has two purposes . Its first function is trivial, 
but required. The number six is typed in column 2 and the rest is 
blank. This indicates to the computer program that there are no more 
input cards. Its second function involves random inputs. This card is 
used to indicate the number of random process (noise) generator cards 
that are to be read before the input process is terminated. Column 2 
has the same information as before. Any pertinent information is con- 
tained in column 5 “ 20 . Columns 21 - 25 must be right- justified and 
they contain the number of random generator cards to be read. 

E. AUXILIARY DECK 

The auxiliary deck is a collection of subroutines and functions 
required by the modules. These routines are general in nature since 
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they can be called by several modules. They are used to calculate such 
things as the properties of air, engine performance, various ratios, and 
to locate values in the many data tables# A brief discussion of the most 
common subroutines is presented. 

1. Boost 

This subroutine calculates the thrust coefficient (CT) and the 
fuel flow rate (FF) during the boost phase. These values are returned 
to module A3 and used to calculate the missile thrust in the body axes 
and the missile weight. 

2. Ramjet 

This subroutine is used during the midcourse, cruise, phase of 
the missile flight. It uses a simplified ramjet model to calculate the 
thrust coefficient and fuel flow rate. This routine is not automatically 
called. The user must designate that he wishes to use it by input ing a 
type 3 card setting FLGRJ(606) equal to one. This then sets up the 
proper stepping in module A3. 

3. Engine 

This is the primary routine during the cruise phase. It is one 
of many routines within the NWC air-breathing propulsion package. This 
entire package is utilized to calculate the engine parameters from inlet 
to exhaust. Again, thrust coefficient and fuel flow rate are eventually 
computed and returned to module A3« The user does not need to supply any 
special input cards to use this subroutine. If no initial value is 
inputed for FLGRJ(606) it is automatically set to zero, which indicates 
this routine is to be used. 

4. Air 

All the properties associated with the air are calculated in 
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this subroutine* This includes the computation of the speed of sound and 
the dynamic pressure. Since the missile does vary in altitude, the 
routine takes this into account as well as the latitude. Once these 
values are calculated, they are returned to module G3<> 

5. Block Data 

This routine contains data tables. These tables cover para- 
meters from aerodynamic coefficients to thermal properties. The total 
data package covers specific bands within the missile operating envelope. 
Data is stored in matrices which includes one, two, and three dimensional 
ones. These data tables are readily available and there are routines 
designed to retrieve this information quickly. 

6. Serch 

This subroutine, along with several functions, is used to 
retrieve information from BLOCK DATA. If the present operating point of 
the missile is not within one of the bands of information, then the 
tables are interpolated. The functions THR2DL, STDLI, STDLIA, and TA3 do 
the interpolation of the tables. Since the tables are of various 
lengths, these functions are very general. 

F. SEQUENCING 

Sequencing is very important in the running of the MCD6DF program. 
Care must be taken to ensure that the modules are processed in the cor- 
rect order at each step. This is essential to eliminate the use of 
obsolete values from the last cycle. An exception to this oroblem is 
the state variables. These are updated simultaneously by the inteerra- 
tion algorithm*- Any module is capable of usin^ the most current value 
of these, no matter what the order of processing. 

To help remedy this problem, module diagrams were devised. Module 



36 



diagrams aid the user in maintaining the proper flow of variables into 
and out of the modules. To design a module diagram start with a box. 

This box will contain all the equations for a specific module. The 
standard procedure for showing inputs and outputs is to use arrows 
pointing in or out from either the left or the right side. The example 
in Figure 5 (Ref* 9) shows this technique. The arrows pointing in from 
the left indicate variable inputs from other modules. Arrows pointing 
out to the right indicate output going to either other modules or 
program output. The last set of arrows points in from the right. These 
show the constants brought in directly or indirectly from the initializa- 
tion module. Since each arrow represents a variable, they must be 
defined. The usual means of labeling the arrows is to use the variable 
fortran symbol and in parenthesis its COMMON ( 3 M 5 ) location. 

Each variable usually has only one COMMON location associated with 
it. In the case of 'state' variables this is not true. State variables 
are defined by four consecutive COMMON locations. The first is for the 
derivative of the variable. The second and the third locations hold the 
lower and upper bounds, respectively, of the integration error. The 
last one contains the variable itself. 

Once all the required module diagrams have been completed, they can 
be combined to get the overall processing order. The usual processing 
order, shown in Figure 6 (Ref. 10) is to start with the geophysical 
group, then proceed to sensors, computers, airframe, and finally 
dynamics. Within each group is a usual processing order and this too is 
shown in Figure 6. This process for determining the program sequence 
will eliminate the use of any obsolete values in the computer run. 
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IV. COMPUTER PROGRAM CHECKOUT 



The basis for the research was the MOD6DF computer program from NWC 
China Lake. The total package received from them consisted of a listing 
of the program and an uninterpr9ted deck of cards. The initial step was 
to input the cards in small groups into the computer and then to examine 
the source listing. This listing revealed that the original deck was 
punched in BCD. This fact was easy to determine since several characters 
were changed (Ref. ll). The library routine NEWD2K was used to translate 
the cards from BCD to EBCDXD. 

Once the translation was completed, an attempt was made to compile 
the new deck. This produced an output which contained many syntax errors. 
These errors were divided into two major groups. Cne effected the use of 
quotation marks in FORMAT and comment statements. The other one, the 
more difficult, effected the DATA statements in the BLOCK DATA subroutine. 

The problem with the DATA statements was due primarily to the differ- 
ence in the compilers used. The compiler at NWC was much newer and 
allowed for the use of more sophisticated inputs. The compiler at NFS 
only allows a data set to start with the first element. This required 
the rewriting of many data groups. To complicate these revisions, a 
limit of nineteen continuation cards is also imposed. These restrictions 
demanded not only the rewriting of many data sets, but also the forma- 
tion of two new ones. 

With the corrections finally completed, the computer would then 
compile the program. The next step was to link all the subroutines 
together. The first attempt was unsuccessful. Inadvertently, the 
subroutine INTR20 had been omitted from the original deck. Using the 
program listing, the contents of INTR20 were typed and included in the 
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main deck. With this addition the program would now compile and link. 

Now that the program would compile and link, attention was turned to 
getting a good simulation. To facilitate this process, two sample 
outputs were obtained from NWC. These outputs cover one missile flight 
which is broken down into a midcourse guidance and a terminal guidance 
s imulation# Hereafter these will be referred to as midcourse baseline 
and terminal baseline, respectively. Using the initial conditions from 
the baseline models, it was hoped that the outputs could be duplicated. 

A. MIDCOURSE GUIDANCE FLIGHT 

The initial run, using the midcourse baseline inputs, revealed an 
overflow problem with the dynamic pressure ( QD( 50 S) ) » Us ing the 
traceback procedures outlined in the Users Manual from the W. R. Church 
Computer Center, the problem was confined to subroutine AIR. The 
problem turned out to be a translation error. The symbol P0 (0 - zero) 
had translated to ?0 in one place and PO in another. This error caused 
the program to use a value left in that memory location from a previous 
run to calculate the dynamic pressure. With this problem remedied, the 
program could progress a little farther. The next stumbling block 
appeared as a divide check. These errors were resolved by introducing 
patches that would bypass a statement that tried to divide by zero. 

Once bypassed, that quantity would be set equal to zero. This was the 
normal procedure of the computer, but it would stop the run after ten 
such errors. Having corrected all these errors, output was obtained 
which covered the desired seventy seconds of flight. 

When the output was examined a major switching problem within sub- 
routine A3 (missile propulsion module) was found. The midcourse 
baseline utilized a simplified ramjet model, but the output was not. 
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ALTITUDE VS. TIME (MIDCOURSE) 



Rather than use RAMJET, the output disclosed that ENGINE had been used. 
The problem was in the logic statement, FLGRJ 0.0. The program was not 
making the desired switch to the simplified ramjet model. To remedy the 
problem, the data card FLGRJ(606) 1 . 0 , was added to the input data deck. 
Another problem was discovered which occured between 14.0 and 14.5 
seconds. During that time period the missile angle of attack ( ALPHA 
(330) ) exceeded ten degrees. When ALPHA exceeds this angle the engine 
is turned off. With the engine off, no thrust is produced causing the 
forward velocity ( VXTP(286) ) to decrease. This limitation was removed 
from the program. After eliminating these problems a flight trajectory, 
Figure 7* was obtained which closely resembled the output of the mid- 
course baseline. A random sampling of the outputed variables were 
compared for exactness. The differences noted were due primarily to 
computer round-off error. 

To further verify the accuracy of the two simulations, three param- 
eter ( TXBA(073)» VXTP(286) ,and YTP(298) )were chosen as representative 
values for the runs. To obtain a feel for the run, these parameters 
are plotted in Figures 8, 9, 10 . Additionally, random time samples are 
tabulated in Tables I, II, III. 

The thrust ( TX3A(073) ) plot, Figure 8, can easily be divided into 
four mission phases. During the separation phase (0.0 - 4.5 seconds) 
the thrust is zero. This is expected since neither of the engines have 
ignited. This is followed by a rapid increase in the thrust. The 
maximum thrust, 27800 lbs, happens during the boost phase ( 4,5 - 9*5 
seconds). At 9*5 seconds the boost motor stops and the transition phase 
occurs for the next 0.3 seconds. During this time the thrust decreases 
rapidly since neither engine is on. Once the engine port covers are 
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clear , the ramjet sustainer ignites. This action initiates the midcourse 
phase (9.8 - 70 seconds) and the thrust remains relatively constant at 
1200 lbs. 

The missile forward velocity ( VXTP(286) ) can be directly related to 
the thrust. During separation no thrust is produced, therefore the 
velocity is zero. As the missile is boosted to Mach two, the velocity 
increases and approaches its maximum value (2800 ft/s) just prior to 
boost engine shut-down. As the sustainer engine port covers clear, the 
velocity decreases slightly. Once the ramjet engine ignites the velocity 
profile is displayed in Figure 9 and it closely resembles that for the 
midcourse baseline. 

The last parameter, y-displacement ( YTP(298) ), was included to 
demonstrate the accuracy of the guidance system Figure 10 shows that 
prior to attaining the cruise altitude, the missile wanders off track. 
Once the guidance system, is activated, it makes the necessary corrections 
to return the missile to the planned flight path. Table III shows that 
the y-displacement corresponds to the baseline case and at the conclusion 
of the run the off-track error is down to 12.03 feet. 

All the information obtained from the new run was checked against 
the midcourse baseline. The results showed that the two simulations are 
within acceptable limits. With this milestone completed, the program 
checkout could proceed to the terminal guidance flight. 

B. TERMINAL GUIDANCE FLIGHT 

To checkout the terminal guidance portion of the M0D6DF program, the 
initial conditions from the terminal baseline were used. This required 
changing about a dozen input cards in the midcourse input data deck. 

The first run attempted turned out successful. This was due to the 
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ALTITUDE VS.' TIME (TERMINAL) 



elimination of all the errors during the midcourse checkout. The result- 
ing flight trajectory, Figure 11, is a very good illustration of the 
desired profile. To check the accuracy of the simulation, the same 
three parameters were chosen. Tables IV and V show that this run is 
almost identical to the terminal baseline. The difference is primarily 
due to computer round-off error. Figures 12 and 13 show the overall 
flight performance of the thrust and forward velocity, respectively. 

For this simulation, the y-displacement ( XT? (298) ) does not 
compare favorably with the terminal baseline. Figure 14 shows that be- 
tween ten and twenty seconds in the flight, XTP develops rapid 
oscillations. By examining Table VI, it is noted that during this 
interval the two simulations differ the most. No reason for this 
discrepancy could be found. These oscillations do effect the flight by 
increasing the time of flight from 25*825 to 29»09 seconds and by 
increasing the miss distance. The imposed time constraints force this 
problem to be overlooked and to direct attention to the accuracy of the 
impact. Since the impact is within 0.5 feet of the target location, it 
is concluded that the terminal guidance simulation works properly. 

With both segments of the MGD6DF computer program working correctly, 
thoughts could now turn towards the actual experimental runs. The 
results of these runs are discussed in the next section. 
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YTP VS. TIME (TERMINAL) 



V. TERMINAL FLIGHT DISTURBANCES 



With the checkout of the M0D6DF computer program completed, many 
ideas were discussed concerning alterations to the simulations. The 
three modifications decided upon were: 

* examine target miss distance when changes were made to the 
initial x-y-z conditions, 

* examine the terminal flight profile, range, and miss distance 
when the cruise altitude was reduced to approximate a sea skimming mode, 

* examine target miss distance when random noise is applied to 
the missile homing seeker. 

To understand the changes and results , one must be familiar with the 
frame of reference. The origin of the reference frame travels from the 
launch platform to the target location at the cruise altitude. This is 
called th9 tangent plane reference system. Displacements in the 
x-direction ( XTP(290) ) are measured from the launch platform in the 
direction of the target. Y-displacements ( YTP(29S) ) are measured left 
or right of the ideal flight path, in the tangent plane. Any vertical 
displacement ( ZTP(306) ) is measured normal to the ideal flight path, 
with down being the positive direction. Using this frame of reference 
the information in Table VII is easier to understand. 

The first modifications demonstrate the effect of changing YTP. 

YTP was increased until a result was reached that was unsatisfactory. 
From the results in Table VIII, the maximum value of YTP was determined 
to be 1800 feet. This conclusion corresponds with the results in 
reference 13 • Since the missile is symmetric, it was also concluded 
that moving YTP either right or left would give the same results. 

The fifth run simulated a drop in the missile's altitude. ZTP was 
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RUN 


XTP 


YT? 


ZTF 


BASELINE 


1735552.5 


0.0 


0.0 


1 


173552.5 


1000.0 


,0.0 


2 


173552.5 


15 C 0.0 


0.0 


3 


173552.5 


1800.0 


0.0 


4 


173552.5 


2000. 0 


0.0 


5 


173552.5 


0.0 


2000.0 











INITIAL CONDITIONS 
TABLE VII 
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RUN 


XTP 


YTP 


ZT? 


IMPACT 

TIME 


BASELINE 


234999.56 


-.0506363 


34999.34 


28.825 


1 


23 ^ 999.81 


.1259258 


34999.836 


29.108 


2 


235000.19 


.11479032 


35000.32 


29.124 


3 


235000.25 


.15263242 


350C0.391 


29.134 


4 


269802.0 


2000.0 


19795.00 


35.0 


5 


234999.63 


.45691 


35000.42c 


30.365 













IMPACT RESULTS 



TABLE VIII 



inputed as 2000,0 feet. As the simulation progressed, the missile 
developed the necessary commands to climb and regain the desired cruise 
altitude (Figure 15) • It then commenced its terminal dive. The target 
was detected, acquired, and impact followed. This initial condition 
modification produced the largest miss distance (0,437 ft.). 

Figure l6 was presented to show the relative locations of the miss 
distances for each change. Table VTII reveals that up until 1300 feet # 
changes to YTP created very little variation in the miss distance. It 
should be noted that the larger the displacement became, the greater the 
time until impact. The increase in time was necessary to allow the 
missile to acquire the target and then compute the required acuator 
commands to impact the target. This concluded the investigation of the 
effects of initial displacement error on target miss distance. 

The idea behind changing the cruise altitude was to try and simulate 
a sea skimming profile. RAMJET contains tables that allow for four 
cruise altitudes. Of these four altitudes, the lowest (500 ft) was 
chosen even though it is high for a sea-skimmer. To run this simulation 
only two input data card changes were required. HQ(4l4) and KREF(50l) 
had to be set equal to the desired altitude. The resultant flight path 
is shown in Figure 17* 

The simulation produced an error-free output, but at first glance 
the results appeared unacceptable. Still trying to compare miss dis- 
tances, the downrange distance (XTP) was found to be 197916 feet. 
Comparing this to the terminal baseline (234999) resulted in an extreme 
error. It was then realized that the missile must expend more fuel at 
this altitude to attain the same speed. Therefore, the results could be 
correct. As further verification of the correctness of the run, it was 
assumed that the missile weight at impact should be fairly equal. The 
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ALTITUDE VS. TIME (TERMINAL - ZTP) 
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ALTITUDE VS. TIME (SEA- SKIMMER) 



simulation impact weight was 1126 . 3 ^ lbs while the terminal baseline 
weighed 11 13*12 lbs. Since no conflicting information could be found 
in reference 12, it was concluded that this simulation was correct. 

The final flight disturbance investigated was the addition of random 
noise to the missile homing seeker. This simulation proved -to be unsuc- 
cessful due to invalid input format. The information describing the 
process is brief and difficult to understand. Further information from 
NWC China Lake is required* to be able to successfully run simulations 
with random noise. 

With the completion of these modifications, many unanswered problems 
and questions still exist. The two most critical problems concern the 
input of random noise generators and the missile exceeding the angle of 
attack limitation when using ENGINE for cruise propuls ion 0 However, the 
simulations did showthat the MCD6DF program runs correctly at alternate 
cruise altitudes and with initial displacement errors. 
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VI. CONCLUSION 



The MOD6DF computer program from NWC China Lake was converted to 
operate on the IBM- 36 O computer at the Naval Postgraduate School, The 
program would only function properly when using the simplified ramjet 
model. When this model was not used, the missile angle of attack 
exceeded the maximum limit of ten degrees. This error caused the ramjet 
engine to flame out. 

Target impact errors for the terminal guidance problem were 
investigated when the initial displacement was modified. These modifi- 
cations demonstrate the missile 1 s accuracy when removed from the ideal 
flight path. The results also point out that if the target falls within 
the seeker search pattern, target impact will inevitably happen. 

Random noise generation is possible with the M0D6DF program. 

However, more information discussing the parameters required to develop 
the noise is necessary. Additionally, sample noise inputs should be 
obtained that reflect the alteration to the noise subroutines. 

This research involved the preliminary investigation of the M0D6DF 
program. The program has many possible areas, concerning guidance and 
control of tactical missiles, which could be developed for future study* 
These areas not only include the unresolved problems encountered during 
this research. Additionally, reference 14 and 15 contain many examples 
of possible advanced guidance concepts. 
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ALCULATE ADDITIVE CRAG CUE TO SUBSONIC SPILLAGE 
TASUE = N/ ( (N*10.2095 1-11 .5 652 ) 

DSUe = 2 . 4 <l.-ETASUBM (l.-ACR) 
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100 CALL EXIT 
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